#########################
data.2000<-dget("data2000")
Yd<-data.2000[[1]]
Xd<-data.2000[[2]]
Xs<-Xr<-data.2000[[3]]

source("gbme.r")
n<-dim(Yd)[1]
Y<-Yd

# call the main estimation package, gbme
gbme(Y=Yd,Xd=Xd,Xs=Xs,Xr=Xr,k=3,fam="binomial",N=matrix(1,n,n),NS=100000,awrite=T,bwrite=T, directed=T) 
